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We demonstrate that graphene-based photonic superlattices provide a versatile platform for elec¬ 
trical and all-optical control of photonic beams with deep-subwavelength accuracy. Specifically, by 
inserting graphene sheets into periodic metallo-dielectric structures one can design optical super¬ 
lattices that posses photonic Dirac points (DPs) at frequencies at which the spatial average of the 
permittivity of the superlattice, e, vanishes. Similar to the well-known zero-h bandgaps, we show 
that these zero-e DPs are highly robust against structural disorder. We also show that, by tuning 
the graphene permittivity via the optical Kerr effect or electrical doping, one can induce a spectral 
variation of the DP exceeding 30 nm, at mid-IR and THz frequencies. The implications of this wide 
tunability for the photonic Zitterbewegung effect in a vicinity of the DP are explored too. 


Introduction. One of the most daunting challenges 
one has to overcome in order to develop ultra-fast ac¬ 
tive photonic nanodevices is the extremely narrow tun- 
ability of physical parameters of optical media. For ex¬ 
ample, while the conductivity of semiconductors can be 
varied by orders of magnitude with the help of impu¬ 
rity doping or thermally, the index of refraction of op¬ 
tical media, which is the main parameter defining their 
optical properties, can hardly be tuned by more than 
0.5%. Several approaches that address this key prob¬ 
lem have been proposed, including the resort to plas- 
monic materials whose optical nonlinearity is enhanced 
by strong optical fields associated with the resonant exci¬ 
tation of surface-plasmon polaritons (SPPs) M , the use 
of strongly anisotropic materials whose optical response 
can be controlled via the polarization of the excitation 
field Hi, reconfigurable metamaterials j§], embedding 
liquid crystals into a photonic-crystal matrix jjj, and tun¬ 
ing the material refractive index via carrier injection [§] . 
While promising if targeted at specific applications, the 
efficiency of these methods is generally limited by high- 
power requirements, slow optical response, and reduced 
spectral tunability. 

A recently introduced photonic platform that contains 
all the ingredients needed to achieve ultra-fast broad op¬ 
tical tunability at the nanoscale is graphene M- Com¬ 
pared to the SPPs in metals, graphene plasmons are 
localized within much smaller regions, provide a much 
larger optical near-field enhancement, and can pass sig¬ 
nificantly longer distances EMI. Equally important, 
the conductivity of graphene and, implicitly, its electri¬ 
cal permittivity can be varied substantially, locally and 
ultra-fast, by means of chemical doping or gate voltage 
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EM3- For example, changing the Fermi level by a 
mere 1 eV causes the relative variation of the real part 
of the permittivity, at the 10 pm wavelength, by 250%. 
Depositing graphene onto the surface of specific mate¬ 
rials, which may be both dielectrics and metals, a vari¬ 
ety of graphene-based layered nanostructures have been 
designed, and their implementations in functional de¬ 
vices have been demonstrated El- The most salient 
feature of these photonic super lattices is that, by in¬ 
corporating graphene into their periodic structure, the 
band-structure and thus their optical response can be 
dynamically tuned. This key property, in conjunction 
with relatively low loss and extremely tight confinement 
of the graphene plasmons, opens up promising applica¬ 
tions of tunable graphene photonic superlattices (GPSs) 
at mid-IR and at THz wavelengths. In addition to these 
important practical considerations, GPSs can also play a 
major role at the fundamental level. Specifically, it has 
been suggested that electron-wave dynamics in graphene 
heterojunctions are in many ways similar to the propaga¬ 
tion of optical beams in photonic superlattices contain¬ 
ing positive- and negative-index material s 1181. intriguing 
phenomena such as the Klein tunneling Il9l 1201 1 an d Zit¬ 
terbewegung (ZB) effect [Til l26j [ being observed in both 
physical systems. 

In this Rapid Communication we demonstrate that, 
by including graphene sheets into the unit cell of metallo- 
dielectric superlattices, one can readily design GPSs that 
posses electrically and optically tunable Dirac points 
(DPs) and feature the ZB of optical beams propagat¬ 
ing near the DPs. A remarkable property of these DPs 
is that they form when the spatial average of the per¬ 
mittivity of the superlattice vanishes, e = 0, an imme¬ 
diate consequence of this fact being their remarkable ro¬ 
bustness against structural disorder. To add specificity 
to these ideas, we consider, as an example of a GPS, 
the one-dimensional (ID) periodic structure depicted in 
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FIG. 1. (a) A schematic of the graphene photonic super¬ 

lattice: the unit cell consists of a graphene-metal-dielectric 
tri-layer configuration stacked along the x direction. The in¬ 
cident TM-polarized beam propagates along the z-axis, its 
non-vanishing field components being E x , E z , and H y . (b) 
The photonic band structure calculated in the zero-averaged- 
permittivity case, featuring the ID Dirac point. The struc¬ 
ture parameters are e g = —14.48, esi = 12.25, e m = —924.05, 
t g = 0.5 nm, tsi = 604 nm, t m = 8nm, and A = 4.2 pm. 


|4c [eV] 




- 0.05 - 0.025 0 0.025 0.05 
kjTt/A] 


0.2 

^ — t»=8nm 


„0.15 

—— ^=20 nm 


& 0.1 
GO 

O 




Q 0.05 

°0 

^ivWTv* 

5 1.5 

~^V*** v ’ 4 

>5 3.5 

J 

4 


kjk„] 



Fig. [Ha). It has a three-layer unit cell, graphene-metal- 
dielectric, the metallic and dielectric layers chosen in this 
work being silver and silicon, respectively. The disper¬ 
sion relation of the structure is easily found by using the 
transfer-matrix method (see Supplemental Material [27} 
for the derivation of this relation), and is given by 
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where k x is the Bloch wave vector, k z the propagation 
wave vector, t Q , tb, and t c are the thickness of the lay¬ 
ers, A = t a + tb + t c is the period of the supercell , 
a = cos K a cos Kb cos k c , and kj = {oj / c) 2 £jPj — k 2 , 
with j = a, b, c, where a, b, and c stand for graphene, 
silicon, and silver, respectively. For convenience, we here 
define Qij = KiTj , Ki = kit z , and Ti = £iU. By fixing 
frequency lo in Eq. (1), dependence k z = k z {k x ) deter¬ 
mines the spatial dispersion relation for the particular 
frequency. 

Tunability and robustness of the Dirac points. A pecu¬ 
liar property of the band structure of such GPSs is that, 
if it possesses two transmission bands, then the bands 
intersect at the center of the Brillouin zone, k x = 0, giv¬ 
ing rise to the photonic DP [see Fig. 1(b)]. For such a 
two-band configuration to occur, the thickness of the di¬ 
electric and metallic layers should be larger than certain 
critical values. As mentioned above, a DP in our pho¬ 
tonic structure exists provided that the real part of the 
average permittivity of the superlattice vanishes, namely, 
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— IHC {Tg + T m + T d ) = 0. 
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DPs are well known in solid-state physics, being at the 
origin of many remarkable properties of recently discov¬ 
ered materials, such as graphene and topological insula¬ 
tors. In this case they are singular points in the electron 


FIG. 2. (a) The dependence of the wavelength of the Dirac 
point on the graphene’s chemical potential and permittivity, 
(b) The photonic band structure calculated for different wave¬ 
lengths. Parameters are the same as in Fig. [T| (c) The den¬ 
sity of states of a 50-period GPS possessing a Dirac point (at 
tm = 8nm) or a Bragg gap (at t m = 20nm). (d) The den¬ 
sity of states of the GPS in (c) that possesses a Dirac point, 
calculated for different disorder levels. 


energy-band structure, and the aforementioned similari¬ 
ties between electron and photon wave dynamics suggest 
that DPs can exist too in the photonic band structure 
of certain periodic structures, one such example being 
the GPS presented in Fig. 1(a). As we demonstrate in 
this study, the graphene component of the present pho¬ 
tonic superlattice allows one to tune the wavelength of 
the photonic DPs by simply varying the chemical po¬ 
tential of graphene, /z c . More specifically, changing p c 
by external means (e.g., gate voltage, chemical doping, 
etc.), one varies the permittivity of graphene and there¬ 
fore one can tune the wavelength, Ad, at which condi¬ 
tion fHe[£(Ar>)] = 0 holds. These ideas are illustrated 
in Fig. 2(a), where we plot the dependence of the wave¬ 
length of the DP on p c . Note that a spectral shift of 
more than 30 nm can be achieved changing p c by 1.5 eV. 
The variation of the corresponding band structure is pre¬ 
sented in Fig. 2(b). 

Similar to zero-h band-gaps, which are formed when 
the spatial average of the refractive index of a Bragg grat¬ 
ing vanishes, DPs in our photonic structure are particu¬ 
larly robust against disorder added to the system. This 
behavior is expected because structural disorder weakly 
affects the spatial average of the permittivity, hence its 
vanishing, which determines the existence of DPs, per¬ 
sists in the presence of disorder as well. To illustrate 
this phenomenon, we present in Figs. 2(c) and 2(d) the 
density of states (DOS) of periodic and disordered GPSs, 
respectively. In the calculations, we used two superlat¬ 
tices, one selected so as to possess a DP (at t m = 8 nm), 
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whereas the other one is a regular superlattice featuring a 
Bragg gap, at t m = 20 nm. The DOS of these two super¬ 
lattices is extracted from dk z /dk x , see the results in Fig. 
2(c). As expected, while the DOS for the superlattice 
featuring a Bragg gap exhibits a conspicuous wide region 
of forbidden states (DOS = 0), the superlattice possess¬ 
ing a DP shows nonvanishing DOS in the whole spec¬ 
trum. The robustness of the DP against the structural 
disorder is confirmed in Fig. 2(d), where we display DOS 
for the disordered superlattice at several different levels 
of the disorder. It is introduced by assuming random 
fluctuations of the metallic-layer thickness around its av¬ 
erage value, namely, the n-th thickness is t 7 ^ = t m o + 5 n , 
where t m o is the average thickness (we choose t m o = 8 nm 
here), and S n is a random perturbation. We assume S n to 
be uniformly distributed in the interval of [—<5, <5], hence 
the disorder level may be characterized by A = S/t m o. 
The DOSs in Fig. 2(d) clearly demonstrates that, even 
when the disorder level is increased to 80%, the DP is 
preserved, appearing as a spectral peak in the DOS, at 
Ad — 4.2 pm. 

Real-time tunability of the band-structure of the GPS 
provides an effective way to control the excitation and 
propagation of optical beams. Figures 3(a) and 3(b) show 
that, when one varies p c from 0.23 eV to 1.5 eV, the band 
structure transforms from a gapless one, which possesses 
a DP, to a gapped structure, this effect being accompa¬ 
nied by a dramatic change in the beam-propagation dy¬ 
namics. At the DP, a normally incident Gaussian beam 
splits into two, as the upper- and lower-band excitations 
are not resolved, due to the fact that the two bands have 
identical slopes near the DP [see Fig. 3(c)]. By contrast, 
when a gap opens the output pattern exhibits four beams, 
due to the difference in the slopes of the upper and lower 
bands near k x = 0 [as per Fig. 3(d)]. Taking losses in 
graphene and metallic layers into account yields similar 
patterns of the output optical field [Figs. 3(e) and 3(f)]; 
however, as expected, the output beams now decay in 
the course of the propagation. In Figs. 3(e) and 3(f) we 
used the scattering time, t = 90 fs, which was measured 
in an ionic-liquid gating configuration |281 ]. the achieved 
Fermi level being as large as 2 eV 0, [3(| ■ 

Optically tunable Dirac points. In addition to electrical 
means, the wavelength at the DPs can be tuned by means 
of the intrinsic optical nonlinearity of graphene, or of the 
dielectric and metallic layers as well. Such nonlinearity- 
induced tunability of the DP is possible in our superlat¬ 
tices because, as explained above, the DPs form when the 
real part of the average permittivity is zero. Thus, non¬ 
linear optical effects that change the refractive index can 
be used to tune the permittivity in or out of the zero- 
average state, thus providing all-optical means to con¬ 
trol the formation of DPs. The optical nonlinearity of 
graphene is expected to be particularly strong, due to its 
unusually large third-order nonlinear susceptibility [3lj] . 
as well as the tightly localized optical field supported by 
graphene sheets. 

Figure 4 illustrates the creation/destruction of a DP 
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FIG. 3. (a), (b) Formation of a bandgap at the Dirac point 
with Ad = 4.2 gm by varying the graphene’s Fermi level. In 
(a) and (c), permittivities of the layers are the same as in 
Fig. [I] In (b) and (d), e g = —195.02. The real part of the 
permittivities in (e) and (f) are the same as those in (c) and 
(d), respectively, whereas in (e) e" g = 1.28%, e'm = 56.37* and 
in (f) e'g = 4.95*, e m = 56.37*. In all cases, the thickness of 
the layers are the same as in Fig. [I] and the input conditions 
select a TM-polarized Gaussian beam with profile E x (x) = 
exp(— x 2 /w 2 ), where w = 2 Ad- 


under the action of the optical Kerr effect, induced in 
graphene and the neighboring silicon layer. At low opti¬ 
cal powers, starting with a superlattice for which e(A) < 
0, the two transmission bands are separated by a gap, as 
seen in Fig. 4(a). When the input power increases, both 
graphene and silicon are affected by the Kerr nonlinear¬ 
ity, their permittivities varying by 5e g > 0 and 8e si > 0, 
respectively. As a result, the average permittivity of the 
structure, e, gradually increases and the corresponding 
band structure is modified as well. At some particular 
value of the optical power, e becomes equal to zero, so 
that, as shown in Fig. 4(b), the photonic gap closes and 
a DP re-emerges. For the particular case illustrated in 
Fig. 4, fesi = 0.07 and Se g = 0.46. The propagation of 
a TM-polarized Gaussian beam, both in the linear and 
nonlinear regimes, is presented in Figs. 4(c) and 4(d), 
respectively. One can see that, when the power of the 
input beam increases, the beam dynamics changes from 
the typical discrete diffraction to ID conical diffraction, 
which is a manifestation of the formation of the DP in 
the latter case. 
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FIG. 4. (a) The linear band structure, calculated for the same 
parameters as in Fig.[l] except e g = —195.02 (e = —0.15). (b) 
The nonlinear band structure, where the Bloch mode of the 
upper band at k x = 0 is used to induce the nonlinear change 
of permittivities, <5esi = 0.07 and 8e g = 0.46. In (c) and (d), 
the input condition is a TM-polarized Gaussian beam with 
E x (x) = A ex p(—x 2 /w 2 ), with w = 2A, and A = 1 V/m in 
(c) or A = 4.13 x 10 7 V/m in (d). The Kerr coefficients are 
n 2 ,si = 1-86 x 10~ 20 m 2 /V 2 and n 2 , ff = 4.09 x 10“ 17 m 2 /V 2 . 

Tunable Zitterbeuiegung effect. Finally, we point out 
that the control of the DP and band-structure of GPSs 
suggest a convenient way to implement, in the classical- 
optics context, a series of intriguing quantum-mechanical 
phenomena. As a relevant example, here we consider the 
ZB effect, which, in the relativistic quantum mechan¬ 
ics, refers to rapid oscillations (trembling) of a moving 
electron around its average position. The ZB of elec¬ 
trons arises because of the interference of particle and 
antiparticle components of the wave function. The elec¬ 
tron ZB oscillations having extremely high frequencies, 
of about 10 21 Hz, therefore this phenomenon is not ac¬ 
cessible to experimental observations. Recently, optical 
analogs of the ZB were studied in ID [22l426j and 2D j2l| 
photonic crystals. Here we demonstrate that GPSs are 
highly suitable physical systems to study this effect. To 
this end, we show in Figs. 5(a) and 5(b) the beam prop¬ 
agation at two different values of the chemical potential, 

= 0.23eV and [i c = 1.5eV, respectively. Note that in 
both cases the GPS has gaps in the wave-vector space, 
the gap at g. c = 1.5 eV being wider. In the latter case, 
the ZB has a smaller oscillation amplitude and a larger 
frequency, which is consistent with the findings reported 
in Ref. [24J . This behavior is more clearly observed in 
Fig. 5(d), where we plot the evolution of the beam cen¬ 
ter, x c = /+“ |E(x)| 2 xdx/f^ |E(x)| 2 dx. 

Still more important, GPSs provide a reliable plat¬ 
form to implement “nonlinear quantum mechanics”, as 
shown in Fig. 5(c), where we examine the ZB dynam¬ 
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FIG. 5. Beam propagation in the linear regime: (a) g c = 
0.23 eV, (b) g . c = 1.5 eV. (c) The propagation in the nonlinear 
regime, with p c = 1.5 eV. (d) Dependence of the center of 
mass of the beam on 2 . The input condition is a TM-polarized 
Gaussian beam with E x (x ) = Aexp(—x 2 /w 2 ), with w = 12A. 
In (c) A = 1.5 x 10' V/m. The thicknesses of the layers 
are t g = 0.5 nm, tsi = 604 nm, and t m = 10 nm, whereas 
n 2 ,si = 1.86 x 10 -20 m 2 /V 2 and n 2 , g = 4.09 x 10 -17 m 2 /V 2 . 


ics affected by the self-focusing nonlinearity, associated 
with the graphene and silicon layers. The figure shows 
that, while the ZB amplitude and frequency remain al¬ 
most unaltered, the self-focusing nonlinearity leads to a 
spatial walk-off of the ZB, viz., monotonic spatial drift 
of the beam’s center in one direction, in the course of 
its propagation. Note that, as the thickness of the lay¬ 
ers of the GPSs can be scaled so that the real part of 
its averaged permittivity remains equal to zero, one can 
easily extend the ZB effect, and the emulation of other 
DP-related quantum phenomena, into the THz and other 
spectral ranges, different from the standard optical one. 

Conclusion. We have proposed graphene-based pho¬ 
tonic superlattices which should provide a highly effective 
and robust approach for the control of photonic beams at 
the nanoscale. By taking advantage of the strong depen¬ 
dence of the graphene permittivity on the chemical po¬ 
tential, it is demonstrated that, due to the peculiar topol¬ 
ogy of their photonic band structure, specially designed 
photonic superlattices, incorporating graphene sheets in¬ 
terlaid with metallic and dielectric layers, feature large 
electrical and optical tunability. Implications of our find¬ 
ings to applications and studies of fundamental effects, 
such as Zitterbewegung, are proposed too. 
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